METHOD AND APPARATUS FOR INTERPRETING HYBRIDIZED 
BIOELECTRONIC DNA MICROARRAY PATTERNS USING SELF-SCALING 
CONVERGENT REVERBERANT DYNAMICS 



Fldd of the Invention 

The invention generally relates to techniques for analyzing biological samples, such as 
DNA or RNA samples, and in particular to techniques for analyzing the output of a hybridized 
biochip to which the sample has been applied. 
Backgropnd of the Invention 
0 A variety of techniques have been developed to analyze DNA, RNA samples or other 

biological samples to identify diseases, mutations, or other conditions present within a patient 
providing the sample. Such techniques may determine, for example, whether the patient has any 
particular disease such as cancer or AIDS, or has a genetic predisposition toward the disease. 
One particularly promising technique for analyzing biological samples uses a DNA 
5 biochip (or microarray) which generates a hybridization pattern representative of the 
characteristics of the DNA within the sample. Briefly, a DNA micioarray includes a rectangular 
array of single stranded DNA fragments. Each element within the array includes millions of 
copies of identical single stranded strips of DNA containing specific sequences of bases. A 
different fragment of DNA may be provided at each different element of the array. In other 
20 words, location (1,1) contains a different single stranded fragment of DNA than location (1,2) 

which also differs from location (1,3) etc. 

A DNA sample to be analyzed is first fragmented into individual single stranded 
sequences with each sequence being tagged with a fluorescent molecule. The fragments are 
applied to the microarray where each fragment bonds only with matching DNA fragments 
25 already embedded on the microarray. Fragments which do not match any of the elements of the 
microarray simply do not bond at any of the sites of the microarray and are discarded. Thus, 
only those microarray locations containing fragments that match fragments within the DNA 
sample will receive the fluorescent molecules. Typically, a fluorescent light source is then 
applied to the microarray to generate a fluorescent image identifying which elements of the 
30 microarray bonded with fragments of the DNA sample and which did not. The image is then 
analyzed to determine which specific DNA fragments were contained within the original sample 
and to determine therefrom whether particular diseases, mutations or other conditions are 
present within the DNA sample. 

For example, a particular element of the microarray may be provided with fragments of 
35 DNA representative of a particular type of cancer. If that element of the array fluoresces under 
fluorescent illumination, then the DNA of the sample contains the DNA sequence representative 

1 



10 



15 



25 



35 



of that particular type of cancer. Hence, a conclusion can be drawn tliat the patient providing the 
sample either already has that particular type of cancer or is perhaps predisposed towards that 
cancer. As can be appreciated, by providing a wide variety of known DNA fragments on the 
micioarray. the resuhing fluorescent image can be analyzed to identiiy a wide range of 
conditions. 

Unfortunately, under conventional techniques, the step of analyzing the fluorescent 
pattern to determine the nature of any conditions characterized by the DNA is expensive, time 
consuming, and somewhat unreliable for all but a few particular conditions or diseases. FIG. 1 
illustrates various conventional techniques for analyzing the flourescent pattern. Some prior art 
systems utilize two or more of the techniques. It should be noted though that many 
combinations of the techniques are not provided in the prior art. Briefly, the fluorescent pattern 
is quantized to yield a dot spectrogram at step 10. Any of four different techniques are 
employed to identify oligonucleotides represented by the dot spectrogram and to then identify 
mutations based upon the oligonucleotides. More specifically, the dot spectrogram may be 
analyzed using a trained neural network recognizer 12. a statistical decision theory recognizer 
14. a fuzzy/expectation method (EM) clustering algorithm 16 or a rule-based inferencing/truth 
table search 18. 

The results are interpreted to yield mutations of interest at step 20. Then the mutations 
are again processed using either a trained neural network recognizer 22, a statistical decision 
theory recognizer 24. a fiizzy/EM clustering algorithm 26 or a rule-based inferencingrtnith table 
search 28. The resuhs are combined at step 30 to yield a diagnosis. 

Finally, disease confirmation is performed at step 32 by "reduction", in other words the 
disease is confinned by pn)babilistic inferencing. 

One major problem with many conventional techniques is that the techniques have poor 
repeatability. Hence, if the same sample is analyzed twice, different results are often obtained. 
Also, the results may vary from lab to lab. Also, skilled technicians are required to operate tiie 
DNA microarrays and to analyze the output resulting in high costs. One reason tiiat 
repeatability is poor is that the signatures within the dot spectrogram that are representative of 
mutations of interest are typically very weak and are immersed in considerable noise. 
Conventional techniques are not particularly effective in extracting mutation signatures from dot 
spectrograms in low signal to noise circumstances. 

Accordingly, it would highly desirable to provide an improved method and apparatus for 
analyzing tite output of tiie DNA microarray to more expedientiy. reliably, and inexpensively 
determine tiie presence of any conditions witiiin tfie patient providing the DNA sample. It is 
particularly desirable to provide a technique that can identify mutation signawres witi,m dot 
spectrograms even in circumstance wherem tiie signal to noise ration is extremely low. It is to 
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these ends that aspects of the invention are generally drawn. 

One analysis technique for achieving the aforementioned advantages is described in co- 
pending U.S. Patent Application 09/253.789 entitled "Method and Apparatus for Interpreting 
Hybridized Biochip Patterns Using Resonant Interactions Employing Quantum Expressor 
Functions", and is incorporated by reference herein. Briefly, the method of the co-pending 
application operates as follows. The method identifies mutations, if any, presem in a biological 
sample from a set of known mutations by analyzing a dot spectrogram representative of 
quantized hybridization activity of oligonucleotides in the biological sample to identify the 
mutations. A-resonance pattern is generated which is representative of resonances between a 
stimulus pattern associated with the set of known mutations and the dot spectrogram. The 
resonance pattern is interpreted to yield a set of confimied mutations by comparing resonances 
found therein with predetemiined resonances expected for the selected set of mutations. In a 
particular example described in the co-pending application, the resonance pattern is generated by 
iteratively processing the dot spectrogram by performing a convergent reverberation to yield a 
1 5 resonance pattern representative of resonances between a predetermined set of selected Quantum 
Expressor Functions and the dot spectrogram until a predetermined degree of convergence is 
achieved between the resonances found in the resonance pattern and resonances expected for the 
set of mutations. TTie resonance pattern is then analyzed to a yield a set of confirmed mutations 
by mappmg the confirmed mutations to known diseases associated with the pre-selected set of 
known mutations to identify diseases, if any, indicated by the DNA sample. A diagnostic 
confirmation is then made by taking the identified diseases and solving in reverse for the 
associated Quantum Expressor Functions and then comparing those Quantum Expressor 
Functions with ones expected for the mutations associated with the identified disease to verify 
correspondence. If no correspondence is found, a new sub-set of known mutations are selected 
25 and the steps are repeated to determine whether any of the new set of mutations are presem in 
the sample. 

By exploiting a resonant interaction, mutation signatures may be identified within a dot 
spectrogram even in cireumstances involving low signal to noise ratios or. in some cases, 
negative signal to noise ratios. By permitting the mutation signawres to be identified is such 
cireumstances. the reliability of dot spectrogram analysis is thereby greatly enhanced. With an 
increase in reliability, costs associated with performing the analysis are decreased, in part, 
because there is less of a require 

ment for skilled technicians. Other advantages of the invention arise as well. 

Although the method of the co-pending application represents a significant advance over 
35 techniques of the prior art, room for further improvement remains. In particular, it would be 
desirable to enhance the method of the co-pending application to achieve a higher degree of 
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repeatability. 

In particular, repeatability is affected by spatio-temporal degradation of hybridization in 
DNA microarrays implementing botit passive and active hybridization. In bioelectronic systems 
implementing passive hybridization, sources affecting repeatability of analysis results include: 

Stochastic variability in diemical kinetics 

Immobilized oligonucleotide damage during fabrication 

Uneven kinetics during thermally facilitated fluidics reaction 

Post hybridization thermal degradation. Currently biochips are "amplification limited". 
TOs is in large part due to losses during high-termperature hybridization downstream. During 
periods when the sample temperature changes from high to low or low to high, extraneous, 
undesirable reactions can occur that consume important reagents and create unwanted and 
interfering chemicals. Rapid transitions ensure that the sample spends a minimum of time at 
undesirable intermediate temperatures, so that the amplified DNA product has optimum fidelity 
and purity. So current methods rely on excessive amplification to compensate for these losses. 

Oligonucleotide entanglement 

Environmental decoherence due to energy and radiation 

Uneven fluidic catalysis 

Unstable fluorescence and chemiluminiscence maricer binding 
Spontaneous emissions 
Partial bindings 

Anti-aliasing during readout and digitization 
Active hybridization is degraded by 

Capacitive coupling between elements of the inunobilizcd matrix 

Partial bindmgs due to current leakage and uneven conductance 

Ultrascale quantum squeeze effects 

Spontaneous emission 

Nonspecific oligonucleotide trapping 

Chaotic relaxation across the array 

Hence, aspects of the present invention are directed, in part, to providing enhanced 
repeatability of biological sample analysis despite these factors. 

Summary of the Invention 

In accordance with a first aspect of the invention, a method is provided for idertifying 
mutations, if any. presem in a biological sample. The method operates to analyze a biochip 
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output pattern generated using the sample to identify the mutations in the sample. In accordance 
with the method the output pattern is tessellated. A stimulus pattern associated with the set of 
known mutations is generated. A resonance pattern is then generated which is representative of 
resonances between the stimulus pattern and the tessellated output patterns. The resonance 
5 pattern is interpreted to yield a set of confinned mutations by comparing resonances found 
therein with predetemiined resonances expected for the selected set of mutations. 

Ir, an exemplary embodiment, the output pattern is a dot spectrogram representative of 
quantized hybridization activity of oligonucleotides in a DNA sample. The stim«h.s pattern is 
generated based upon Quantum Expressor Functions. The dot spectrogram is tessellated to 
10 match morphological characteristics of the Quantum Expressor Functions and local parametncs 
are extracted. The tessellated dot spectrogram and the stimulus pattern are transformed to a 
metrically transitive random field via phase shlftmg. The resonance pattern is generated by 
iteratively processing the tessellated dot spectrogram by performing a convergent reverberation 
to yield a resonance pattern representative of «sonances between the Quantum Expressor 
15 Functions and the tessellated dot spectrogram until a predetemiined degree of convergence .s 
achieved between the resonances found in the resonance pattern and resonances expected for the 
set of mutations. THe convergent reverberation includes the step of performing a convergent 
reverberant dynamics resonance analysis of the tessellated dot spectrogram using the resonance 
stimulus pattern to identify mutations represented by the tessellated dot spectrogram. The 
20 convergent reverberation also includes the step of performing a convergent reverberant 
dynamics resonance analysis of the mutations using the resonance stimulus pattern to ,dent.fy 
diagnostic conditions represented by the mutations. 

Also in the exemplary embodiment, the convergent reverberant dynamics resonance 
analyses are performed by determining resonance dynamics relaxation values based upon the 
25 tessellated dot spectrogram and tiie resonance stimulus; filtering die dynamics relaxation values 
using ensemble boundary and complete spatial randomness (CSR) filters to yield a second set of 
values; applying bulk property estimators to the dynamics relaxation values «, yield a third set 
of values; evaluating the second and third sets of values to determine a degree of resonance 
convergence; and then determining from the degree of resonance convergence whether a 
30 paralysis of dynamics has occurred and. if so. repeating tiie aforementioned steps. 

In the exemplary embodiments, by tessellating the dot spectrogram to match 
morphological characteristics of the Quantum Expressor Functions and by exploiting a resonant 
interaction employing a resonance convergence check which uses extracted tessellation 
parametric, mutation signattires may be identified within a dot spectrogram with a h.gh degree 
35 of repeatability. By achieving a high degree of repeatability, the reliability of dot spectrogram 
analvsis is thereby greatly enhanced. With an increase in reliability, costs associated wtth 



In accoraan spectrogram representative of 

thet^y The me*od compnsing s.q>s of d«=n.ining r«..»c dyn™« 

rpi^ and electron c microdevices. Principles oi 

Lion analysis at ultra-low concentrations using ultra-high dens.ty pass.ve and/or 

andprotein-DNAinteractionsforn.edicaldiagnosisatthemolecularlevel. 
Apparatus embodiments are also provided. 

Rrij'f DescriDtion «f the Drawings 
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apparent from the detailed description set forth below when taken in conjunction with the 
drawings in which like reference characters identify correspondingly throughout and wherein: 

FIG. 1 is a flow chart illustrating exemplary techniques for analyzing dot spectrograms, 
some combinations of which are in the prior art. 
5 FIG. 2 is a flow chart illustrating an exemplary method for analyzing the output of a 

hybridized DNA micioarray biochip in accordance with the method of the co-pending 
application. 

FIG. 3 is a flow chart illustrating exemplary enhancements to the method of the co- 
pending application provided in accordance with the present invention. 
10 FIG. 4 is a flow chart illustrating an exemplary method for analyzing the output of a 

hybridized DNA microairay in accordance with the method of the present invention. 

Detailed Description of Exemplary Embodiments 

With reference to the remaining figures, exemplary embodiments of the method of the 
15 co-pending application and of the present invention will now be described. The exemplary 
methods will be described primarily with respect to the analysis of mutation signatures within 
output patterns of DNA biochip microarrays, but principles of the invention to the analysis of a 
wide variety of other patterns as well. 

20 Overview of Co-Pendintf Application 

Briefly, the exemplary method of the co-pending application exploits, among other 
features: (a) a novel representation, interpretation and mathematical model for the immobilized 
oligonucleotide hybridization patterns, represented via a dot spectrogram; (b) a new "active" 
biomolecular target detection and discrimination method based on quantum resonance 
25 interferometry, and (c) a new spatial hashing function that yields accurate diagnostic assessment. 

To this end the exemplary method of the co-pending application exploits a 
fundamentally different computational paradigm for mutation expression detection in pre- 
enhanced dot spectrogram realizations. The method is an innovative modification to 
dynamically anayed quantum stochastic resonance (QSR) for discrete system analysis. The 
30 arraying strategy is a fimction of the expression pathway of interest. The method depends on the 
molecular diagnostic spectrum being addressed. Banks of coupled quantum resonators are 
algorithmically designed to significantly enhance signal-to-noise (SNR) performance and fuse 
multiple synthetic renormalized dot spectrogram realizations to better detect prespecified 
biomolecular expression patterns. 
35 Moreover, the exemplary method of the co-pending application exploits an enhancement 

in previous extensions to classical stochastic resonance (SR) and array enhanced SR (AESR) in 
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signal processing and sensor data analysis. Stochastic resonance is a phenomenon wherein the 
response to a sensor, modeled in terms of a bistable nonlinear dynamical system, is enhanced by 
applying a random noise element and a periodic sinusoidal forcing function. SR occurs when the 
SNR passes through a maximum as the noise level is increased. 

Tlius as important aspect of the exemplary method of the co-pending application 
involves the coupling of transformed and preconditioned discrete microarray outputs to a 
mathematical model for a quantum-mechanical dynamical system with specific properties. 
When driven in a particular manner, the coupled system exhibits a nonlinear response that 
corresponds to detection of phenomena of interest. TTie method exploits modulation of 
observables from a "base" (canonical continuous dynamical system), so that a selected set of 
spectral properties match a similar selected spectral properties of a discrete spatial tessellation 
substructure from an amplitude spectrogram derived from bioelectronic observables. The 
method further exploits the concept of convolving a discrete spatial system (derived from base 
mutants of interest) with a continuous asymmetric temporal system to derive a spatiotemporal 
15 input to further convolve with another discrete spatial projection (of an inherently partially 
stabilized spatiotemporal system). 

Hence key components of the exemplary biomolecular detection method of the co- 
pending are: (i) selection of a basis system; (ii) generation of designer Quantum Expressor 
Function (QEF) for coupling with the substrate to be analyzed; (iii) generation of a HamUtonian 
20 to describe relaxation dynamics of the coupled system; (iv) modulation of resonance parameters 
to enforce early resonance; (v) and exploitation of resonance suppressors to verify detection. 

Referring to FIG. 2. initially at step 100. a set of mutations of interest are selected. The 
mutations, for example, may be mutations relevant to cancer. AIDS, or other diseases or 
conditions. At step 101. preconditioner transforms are generated based upon the selected set of 
25 mutations. TTie preconditioner transforms are provided to convert mutation nucleotide 
sequences into expected amplitude patterns in the prespecified microarray representation, given 
a particular biochip layout. At step 102, Quantum Expressor Functions are generated based 
upon the Hamiltonian of a pre-selected basis system. The Quantum Expressor Functions are 
designed to couple the Hamiltonian for the selected basis system to a predetemiined DNA 
30 microarray configuration to permit a resonance interaction involving the output of the DNA 
microattay. Resonance stimulus is generated, at step 106, using the Quantum Expressor 
functions, 

What has been summarized thus far are preliminary steps performed off-line for setting 
up the Quantum Expressor Functions and the corresponding resonance stimulus. These steps 
35 need be perfomied only once for a given set of mutations and for a given DNA microarray 
configuration. Thereafter, any number of output patterns from the DNA microarray may be 
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processed using the Quantum Expressor Functions to Identify whether any of the mutations of 
L p«-selected set of mutations are found therein. Preferably. Quantum Expressor Funcfons 
are pre-generated for a large set of mutations and for a large set of DNA ™croarray patterns 
such that, for each new DNA microarray output pattern from each new patient sample, the 
5 p^sence of any of the mutations can be quickly identified using the predetermined set of 
Quantum Expressor Functions. In general, the aforementioned steps need be repeated only to 
update the Quantum Expressor Functions to accommodate new and different DNA microarray 
patterns or to if new mutations of interest need to be considered. 

At step 106. an output pattern (referred to herein as a Dot Spectrogram) is generated 
10 using a DNA microarray for which Quantum Expressor Functions have already been gene^^d^ 
At step 108, the dot spectrogram is preconditioned to yield a dot spectrogram tesselaUon (DST) 
to permit exploitation of a resonance between the dot spectrogram and the Quantum Expressor 
Functions The actual resonant interaction, which involves convergent reverberations, .s 
performed at step 110 until a predetermined degree of convergence is achieved. Once 
15 convergence is achieved, a resulting resonance pattern is processed at step 112 to ident.fy any 
mutations represented thereby.As will be described below.steplll is Hindered trivia.^ 

of the aforementioned resonant interacUon which is based upon Quantum Expressor Function 
already conflated with the pre-selected mutations. Hence, no complicated analysis is ,.qu.red 
to interpret the resonance pattern to idendfy the mutations. Next, at stepll4. the mutat.ons.^ 

20 mapped to cotresponding diseases and conditions to thereby identify any diseases or cond.t.ons 
that the patient providing the sample being analyzed is afflicted with. Again, th.s .s a fa.ly 
trivialstep. Finally,atstepll6,diagnosticconfirmation is preformed to verifythatthed.seases 
or conditions are present in the sample. This is achieved by starting with the found diseases or 
conditions and then performing the steps of the method in reverse. 

25 Each of the aforementioned steps are described in detail in the co-pending apphcation 

and the detailed description thereof is not repeated herein. 

n i.w nf the Ex '-^pi'^ry Method of th.^ Present Invention 

The present invention is directed, in part, to improving the repeatability of the method o 
30 the co-pending application by tessellating the dot spectrogram so as to match morphologtcal 
cha^cteristicsofthe Quantum Expressor Functions andby using extracted local ^^^^ 
part ofaresonance convergence check during the resonance interaaion.Tl.eseaddm^^^^^ 
have the advantage of establishing uncertainty bounds which permit method repeatabthty to be 
enhanced and quantified. 

3, FIG. 3 illuara^s «*«>c«n»« <o .=l"iq- of ™. 2 P^vded by .h. preen. 

.tag s«ps of*. »ch»K^ of FIG. 2. m s^s omo. 1 .ppcamg 
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in FIG 3 may be the same as those of FIG. 2 and will not be redescribed. Like reference 
numerals, incremented by one hundred, are employed to represent the repeated steps. 
Briefly the enhancements of the present invention entail five major steps: 
. Preconditioning the hybridized array output pattern (i.e. the dot spectrogram) by 
5 afuzrytessellationandcouplingthepreconditionedoutputpattemwithacanonical^stcmw.th 

aftereffect and memory properties (step 218); 

Estimation of amplitude wanderings or dispersion (step 220); 
. Implementation of resonance interaction by integrating partial or subthreshold 
resonances using phased array enhancement operator resonance dynamics (step 224) w.th 
,0 additional resonances synthetically induced to accommodate the possibility for the presence of 
single-point and two-poim mutations around the mutation-centered pixels; and 
Resonance recombination (step 226). 
A combination of one or more these enhancements are superimposed on the techniques 
described in the co-pending patent application to address specific sources of hybridation 
15 degradation, device imperfections and protocol variability in the analysis p«cess to thereby 
enhance repeatability. 

Referring now to HG. 4. initially at step 300. a set of mutations of interest are selected 
and at step 301 preconditioner transforms are generated based upon the selected set of 
mutations. At step 302, Quantum Expressor Functions are generated based upon the 
20 Hamlltonian of a pre-selected basis system. Phase shifted resonance stimulus .s generated, a 
step 304, using the Quantum Expressor functions. Grouping stimulus is also generated, at step 
305 

steps 300 - 305 are preferably performed off-line to set up the Quantum Expressor 
Functions and the corresponding resonance stimulus and grouping stimulus and need not be 
25 repeated other than to update the Quantum Expressor Functions to accommodate new and 
different DNA microarray patterns or if new mutations of interest need to be considered^ 

. At step 306. a dot spectrogram is generated using a DNA microarray for which Quantum 
Expressor Functions have already been generated. At step 307. the dot spectrogram .s 
tessellated to match morphological characteristics ofthe Quantum Expressor Functions y.eW^^^ 
30 adotspectrogramtesselation(DSTO. At step 308. local parametrics of the tessellated .mage are 
extracted. Then, at step 310, an amount of amplitude wandering is determined and compared 
with pre-determined allowable generator function limits. If. at step 310. the amphtude 
wandering is not within the allowable generator function limhs. then execution p^^^^^^^ 

312 where the tessellated dot spectrogram to match spectral characteristics of the Quanmm 
35 Expressor Functions. Steps 308 and 310 are repeated untU the amplitude wandering is found to 
be within the pre-determined limits at which point execution proceeds to block 314 wherem a 



resonance interaction is performed between the tessellated, renormalized dot spectrogram and 
the phase shifted resonance stimulus generated at step 304 and the group stimulus generated at 
step 305 to identify any mutations represented by the dot spectrogram. 

•nie actual resonant interaction, which involves convergent reverberations, includes the 
5 following sub-steps also shown in FIG. 4. At step 316, a resonance dynamics iteration is 
initiated which includes the use of ensemble boundao' and CSR operators (step 318) and the use 
of bulk property estimators (step 320). The ensemble boundary filters. CSR filters and bulk 
property estimators are applied to the tessellated, re-normalized dot spectrogram in combmation 
with the resonance and group stimulus. The resulting filtered dot spectrogram is then evaluated 
10 to determine a degree of resonance convergence to one or more of the set of predetermmed 
mutations, at step 322. The degree of convergence is evaluated, at step 324. against a Lindbald 
condition and, if the Lindbald condition is not met, the system is deemed to be subject to 
paralysis of dynamics and execution proceeds to step 326 wherem possible hfacel death .s 
compensated for by increasing a time scale for the iteration Initiated at step 316 and then 

15 repeating the iteration. 

Here it should be noted that mutation death and paralysis of dynamics are different 
concepts. THe mutation death check is a conditional check . If this check shows that a resonance 
is not possible for a specific mutation resonance centered (MRQ-hixel then the iteration .s 
terminated and block 314 is exited. But failure of resonance dynamics is not sufficient to 
20 conclude that a specific mutation is absent Indeed, if the "hixel death" check fails, that imphes 
that resonance could be still obtained In a downstream iteration. 

If. at step 324, the iteration has converged and no paralysis has occurred, then one or 
more muUtions have been reliably identified. Execution proceeds to block 325 wherein another 
resonant interaction is performed to identify particular diseases represented by the mutations^ 
25 Briefly, at step 326. a resonance dynamics iteration is initiated which includes the use of 
ensemble boundary and CSR operators (step 328) and the use of bulk property esthnators (step 
330) The resulting filtered dot spectrogram is then evaluated to determine a degree of resonance 
convergence, at step 332. The degree of convergence is evaluated, at step 334, agamst a 
Lindbald condition and, if the Lindbald condition is not met. the system is deemed to be subject 
30 to paralysis of dynamics and execution proceeds to step 336 wherein a time scale for the 
iteration initiated at step 326 is increased before the iteration is repeated. If the Lmdbald 
condition is met, then diseases corresponding to the mutations found using step 314 have been 

reliably identified. r^r^ 
Processing in accordance with step 335 depends on the biochip. The flowchart of HG. 4 
35 illustrates a general form of the method requiring precessing during step 335. In other 
implementations, this step is trivial or can be eliminated entirely. In this regard, the overall 
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method is implemented at two levels of abstraction, depending on how well the disease 
genomics is understood. Detection of a specific mutation is necessary and sufficient to conclude 
expression of a specific gene. But the expressed gene may not be the one to conclusively 
identify the disease. Then another level of abstraction is invoked wherein the method is applied 
inferentially by expanding the gene expression circuit or gene expression tree to determine if 
there is evidence that all expressed genes eventually lead to one that conclusively identifies the 
disease. So step 335 operates on the results of step 314, such that all identified mutations are 
used as an input to determine if the complete expression pathway for leading up to the point that 

a disease can be concluded. 

If the biochip is so designed that the mutations corresponding to al! intermediate 
expressed products, from any disease starting point can be captured by resonance output of step 
314. then sub-steps within step 335 and subsequent steps 340 and 342 can be circumvented. If 
not, clustering step 340 and geometric hashing step 342 are provided to identify that an 
expression pathway is present that trivializes the disease conclusion Kstep 344). 

Once the diseases are identified, clustering properties are evaluated at step 340 to 
selectively eliminate oligonucleotides representing possible diagnoses based on morphological 
filtering of subthreshold resonances and any subsequent reccntering (i.c..the inverse of 
dispersion). Steps 314 - 335 produce a cluster of sub-threshold resonances. Step 340 is a 
reverification such that all induced resonances are present in the target sample and not a 
20 manifestation of multiple rescalings and synthetic SNR enhancements. 

Then a hashing projector is applied at step 342 to order the mutations. A diagnostic 
decision is then rendered at step 344 by examining the order of the mutations and comparing the 
mutations with a table identifying corresponding diseases. 

Thus, the output of block 314 represents all hixcls that identify complementary 
25 oligonucleotide bindings in the biological sample being analyzed and this represent "mutations". 
TTie output of block 335 comprises a set of expressed genes that are associated with a particular 
pathogenic pathway and thus represent a preliminary "diagnosis". Further analysis of the 
pathogenic pathway provides a set of possible diseases, if any. This decomposition is motivated 
by scaling the computation to answer three questions: 

Does the currem condition imply predisposition to a specific disease? 
What is the likelihood of presence of a disease? 
What is the set of all possible diseases that may be concluded from the target 
sample, given a specific genomic encoding implemented by the biochip? 
In any case, if the diagnostic decision rendered at step 344 is affirmative, then the 
35 diagnosis is output If the diagnosis is negative then, at step 346, a determination is made as to 
whether there are any alternative mutations, not within the initial set of mutations selected at 
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step 300, that could be present within the sample. This detemination is made by examining a 
table listing all possible mutations. If there are alternative mutations, then the process is 
repeated from step 300. If not, then a signal is simply output indicating that no mutations were 
found in the sample. 

5 Now details of the steps of the new method will be provided. Details regard steps 

already described in the co-pending application will not be repeated herein. 
Mutations Sets 

TTie mutation set of Interest generated at step 300 is selected by identifying 
oligonucleotides representative of the {Z) mutations of interest. Each oligonucleotide is 
10 represented by Mm) which is given by [«oa....aJ. where a = {A,C,T,G} base associated with 
each array cell [a,b] where 10<*<25. The entire set of unique oligonucleotides denoting 
mutations of interest. A(l,m), is given by [5o5,...6J where 5 = {A,C,T,G} length|Al = 
length 1 v I , and 0<1| A - M k. and the designed in vKLm) oligonucleotide sequence is a 
perfect complement to only A(l,m) for all l.m. 
15 As part of step 300, an oligonucleotide table is generated which contains the 

oligonucleotide sequences associated with each mutation of interest identified by row and 
column location (ij). The oligonucleotide table is provided for subsequent use at step 312 to 
map locations within the dot spectrogram wherein resonance occurs at step 310 to 
oligonucleotides such that mutations present in a sample being analyzed are easily identified. 
20 Also as part of step 300, a mutation table is generated which contains the diseases associated 
with each mutation of interest The mutation table is provided for subsequent use at step 314 to 
map mutations identified at step 312 to specific diseases or other medical conditions such that 
the diseases can be easily identified. 

25 Basis System for Quantum Ey pressor Functions 

The selection of the basis system and the generation of the QEFs based thereon 
depends, in part, and the characterisitcs of the DNA microarray. In the exemplary embodiment, 
the DNA microarray is an N by M DNA chip array wherein an element of the array is referred to 

herein as an "oxel": o(ij). 
30 The pre-hybridization microarray (PEBC) is expressed as: 

PEBC = XZoOW). where N and M refer to the linear (row and column) 
dimensions of the 2-D microarray. 



The numeric value associated with each oxel is given by: 
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where [al=[A | C 1 T | G] take the values [0 1 1 1 2 1 3 ] respectively. 
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An element of the dot spectrogram is referred to herein as a Wxe/; h(ij). 
A spin boson basis system is selected for use with this type of array. Other basis system 
be appropriate for either the same or other microarray configurations.. 



Quantum Px pressor Functions 

n,e QEF is generated at step 302 based upon the spin Boson basis system by first 
,0 calculating the Hamiltonian for the system, calculating hamionic amplitudes |P^ for the 
Hamiltonian. generating an order function (OF), measuring entraimnent entrainment states of 
the OF of the ground truth and finally modulating the OF of ground tmth to yield the QEF. 

The QEF's generated at step 302 are converted to a phase-space representation. Also, .f 
the output of the hybridization chip is not in phase space then it is converted as well. The 
,5 conversion is performed using phase embedding operator. T, described in the co-pendmg 
application. Results associated with combinatorial Hopf Algebra are used to contain amphtude 
dispersion due to loss of hybridization. A special case of quantum random walk. Gelfand- 
Naimark Segal (GNS) construction is used to disperse group stimulus. Note that coproduct 
construct of the Hopf algebra plays the role of "sharing ouf ' possible explanations of a fact. Tl^e 
20 GNS dispersion of QEF is implemented using an approximation: 

<p^^iampyector) = I/, flt/;' where U,=e-"'' 

where H denotes the hamiltonian for the coupling spin boson system. 

25 

n^^ration andT'-^^^lation of the Dot Spectro.^ 

As noted a dot spectrogram is generated at step 306 for a sample ficm an N by M DNA 
chip array wherein an element of the array is an o(iJ). A 6-a manufacturing process 

accuracy in microarray design is assumed. Each array cell amplitude is given by 0(i.j) for .: 1 
30 to N, and j: 1 to M. Let M/(ij) denote the a priori known oligonucleotide given by [cxoa,...akl. 
where a = {A.C,T.G} base associated with each array cell [a.bl where 10 <^:^ 25. The 
complimentary strand, derived from unknown sample is denoted by ^HJ) • 

The post-hybridization microarray is treated mathematically using the machmery of 
equations with aftereffect. Each hixel given by «l>(ij) is represented as a cluster of dynamical 
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i»oo.«.,ybc™d.Th»s(CB)„w, ^["Bl«, *tPBl». -T«,.«*».0.00«1^ 

do. sp.«n>g™. Wo) is *« «s»U««i » dettrmin. idolized en3»,bte 

boundaries for forcing dowmneamtesommt action. 

TypWiy. in signa, processing appilc^ions, hi^ pass or pas, s^».l »l«nng » 

i„p,e.eL .0 enhance SNR in DS m«rix. Alt»,«e n»«^ apP>y • comb,^ of 

:uLoro.Kered,eae»c,ionfi,«,sa„..oenha^»'««>'-''-'*"7*^^ 
^ ridi»,ionco»ce„>r..lonfto.-»» of*. ^i«en,c.,b. These SNRenhancen-en.^^^ 

lever on,, wi* posUi» or Since S»« i. gene.. 

«„,,,.-,owurge.DNAco.^«i»,*-'-^'"''f'«-^'*~-""'^"" 

""11llspe,,br™edbyperfo,n.^..aien.refoc.in.a.drescaiin..3desc^^ 
*e co^nding app,ic«ion. .n .».ma,iv., a Dirichie, .esse„a«o. operator or a Deiaona, 
triangulatioo operaor are applied to t«»ella«i *» do. spectragran^ 

U. tessellated image is treated as7;«So:il7 -do," 8««- « l-^^ 
^i.,ed»i.has,ngolar(de.e™inistic,,homogeneo»s,i.,s.ati.na,y)f.eU»e*.n.e4 

, The par^netric of most interest is the integrated density of states, given by 



Where = VX-;^ Sin' ^ 
25 whence 



30 



wheZ is the number of eigenvalues to the system (random field approximation). 

Thisiscomputedforeachtessellationregioninthedotspectrogram _ 
Amplitude wandering is determined using Pahn generators as descr.bed m the co 

. \ln m Palm generators exploits the notion of generator functions to capture 
pendmg patent #1. The Palm gen ^^^^^.^^ ^^^^^.^ 

stochastic variability in hybridization bmding efficacy. The exemplary 

dmws upon n^lts in stochastic integral geometry and geometric probabU.ty theory. 



"Amplitude wandering estimate" that bounds the hixel amplitude dispersion due to total 
hybridization losses, is computed using Palm generators over the globally re-scaled dot 
spectrogram to capture amplitude wanderings and transitions at element, neighboring pa.r and 
local ensemble levels. Step 310 provides a measure for each mutation-r«:ognizer cente^d 
5 (MROhixel thatisinvarianttolocaldegradation.Themeasureisexpressedviatheform 

m(Z)= [/(r)cfe 

where Z denotes the set of mutations of interest In other words, we determine the 
Action f(z) under the condition that m(z) should be invariant with respect to all dispersions ^. 
10 Also up to a constant factor, this measun. is the only one which is invariant under a group of 
n^otions in a plane. In principle, we derive deterministic analytical transformations on each 
MRC-hixel., that map error-elliptic dispersion bound defined on tR^ (the two dimension 
Euclidean space - i.e., oxel layout) onto measures defmed on S«. The dispersion bound is given 
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by the form 



Recall that Palm distribution, U of a translation W invariant, finite intensity point 
process in <R" is defmed to the conditional distribmion of the process. It is expressed in terms of 
20 a Lebesgue factorization: 



EpN*= ALn X n 



Where n and A completely and uniquely determine the source distribution P of the 

25 translation invariant point process. ^ term E. denotes the first moment measu« of the 
point process and L^isaprobability measure, in patent^l we described how to computenand 

A which can uniquely encode the dispersion and amplitude wandering associated with the MRC- 



hixel. 
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,„ m mv«>tion w. r.ta the ««»g «s«n,ptto Itet P.ta g.n«.««. n «>d A, 
a„ s^.. of ^ochasUclW in do. spec«^ o«p* Si... h>*ridi^i«n - » 
untax-™ »d unp^dicabl. ways, «. n.«i » modifr «« « P'^>>^ f^-'"' 

^W^. in od« worfs .he 8»««n, arc <o ^'t^ as oppc=«. » » pom. 



function. 
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(n a n m) specifies a continuous probability density function for amplitude 
wandering in the m-th MRC-hixel of interest where the terms denote: 
oligonucleotide density per oxel p,,,,. PGR amplification protocol (o^, fluorescence binding 
5 efficiency (^™) and imaging performance (x.„). Previously we requited a preset binding 
dispersion limit to be apriori provided to compute A, given by the second moment to the function 
atSNR = 0. 

Nondispersive n is computed using U=@*P where 

10 And T, and t, represent normalized hybridization dispersion limits (typically preset to 0.1 

and 0.7 respectively to assume losses between 10% - 70% hybridization. 

Preconditioned dot spectrogram is represented by where function 

l/(l+exp(p(...))) was used to express the underlying known and stationary point process. 

The latter assumption is relaxed in this method and determination of whether the 
nplitude wandering is within allowable generator fimction limits is achieved by: 



amc 
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rViML < n. < 



where 

no denotes the nondispersive generator. 
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^, ^ provide the Laplace characteristic functional of the Poisson random field 
associated with each source of hybridization degradation. The contributions are estimated using 



^,=c,.(det5.)->^A,'/^ 
30 where c is gain constant and a denotes a nonrandom matrix 
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solution of the following variational problems: 



TT.e differential operator for the metrically transitive field for convoling the uncertain^ 
parametersisdcnotedbyAo. Population and solution ofthe above equation requires eshntates 
fortheforwardsensitivitymatrixofvariablesimpactinghybridizationdeg,.dat.on. 

aenonnal^tion at step 312. if necessary, is performed on the tesseUated .mage . 
^rther match spectral properties of the sHmulus pattern 1.e re— .n of 
spectrogram is achieved by rescaling the dot spectrogram in .he .nterva. -l" 
rJon proceeds in the phase space which is why we transformed the system to the 
metrically transitive random field. 

p^nn.nt Interactic" ^" '^^"tifv Mutations 

between d.e QEP and the tessellated. 
25 ncrmalizeddotspectrogramisperformeduntilapre-selecteddegreeof.^^^^^^ 
Resonance dynamics relaxation values are calculated at step 316 as follows. 

A closed-form convolutionless evolution equation is given by: 

30 V (t) =0it,r) i^J) 
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where 



18 



whe« both depend upon the nonnalized DST, (i.e.. initial state) at time x - post- 
hybridization but pre-conditioned state. And 

5 

^ (t,T) is ihe unitary evoluter e 

DST 

for the hamiltonian specifying the after affect basis system dynamics. 

10 Also, 

^(/.r) = limi[<^(/ + ^.T)-l] 

In practice a small e .typically lO"* is used. 

SO if theoretical convergence time is r. (outer convergence cycle timOandchoosingt: 
T+ To, then: 
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Dsr, ' 
20 and 



25 and Jt is a time-ordering operator. 

The dynamics relaxation values are then filtered at step 318 using ensemble boundary 
and CSR filters (higher order Poisson kernel) as follows: 
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<p {fye''"^ where 

DST, 
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_ where rSO. 

l-2rCos+r 



The bulk property estimators of step 320 are applied to the dynamics relaxation values as 
follows: 

where t is the discretization step. 

•n,e above expression provides an estimate of when a geomteric motion embodied by the 
convolutionless equation, is no longer a plausible resonance candidate. 

TOs is the closed form for an expression at which the coupling between DST and the 
niicroarray is broken and a coupling with a nonlinear information filter (NIF) is established. In 
essence, the system forgets any initial correlation and tends to a Lindbald condition. 

The resonance convergence is determined at step 322 as follows: 



U(t+1)-U 



Log 



[/(/)- L/| 



Si 



The system oscillates if no convergence is reached. If increasing the timescale x-times 
20 (~5) does not meet the condition, then the mutation is deemed to be absent. 

It should be noted that, unlike the technique of the co-pending application, in the present 
invention the absence of resonance over a maximum interation count does not imply absence of 
resonance, m reason is that both the dot spectrogram and the QEF are dispersed, i.e.. the SNR 
is reduced over an individual hixel, but is in fact increased over an ensemble. So the 
25 convergence decisions are made by cascading the inner loop reverberations as opposed to a 
single reverberation. So two timecycles are used for the convetgence analysis: 

(a) time cycle over which hyperfme resonances are tracked, detected and used as a dec«.on 
mechanisms to continue or stop the interation; 

(b) time cycle over which the absense of mutation is actually concluded. This is done by 
30 implementing a local maxima over output of previous step and then reintegrating.Tl.e 
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method essentially accumulates partial resonances and then applies the same resonance 
equation to the rescaled and renormalized partial stage. 
This process can be analytically be represented as: 



L 



C(T,) 



DST OUT, 



> 0 



0 



otherwise 



where c. c, and a,* thresholding constants that are used to detect subthreshold 
resonances. Also. c,»C2 »c,» 1 /[amplitude resolution]. 

Also X, and refer to the inner and outer integration timescales. In an implementation 
they refer to the iteration conter at which the integration loop is terminated, exceeded or exited. 
,0 Typically termination counter is set to one thousand steps with timescale of the order of ten 
nanoseconds for inner step and microsecond for outer step. So effective device convergence 
time is within onen hundred milliseconds for the entii* computation. 

In this regard, if the Lindbald condition is not achieved and verified the dynamics .s 

considered paralyzed. 
15 TTie paralysis of dynamics actually implies that coupling between 

(a) the Quantum Expressor Function, 

(b) the tessellated and normalized dot spectrogram , and 

(c) the convolutionless carrier 

is too weak to exhibit a nonlinear resonance. The physical interpretation is that the coupled 
20 system exhibits "frustrated dynamics" which enhances and impedes resonance reaction at the 
same time. So the actual output takes the form of white noise over several hixels which 

oscillates. . 
The detection of oscillation occurs when the spectral radius for the convergence crrtena 

oscillates between limits [e„ e. ] and does not tend towards O.Tltis may be verified by trackmg 

25 the spectral radius zero crossing with respect to the lower bound e.. If the zero crossing 

frequency exceeds a present number (e.g.. 10) in this implementation, the dynamics is deemed 

paralyzed. , , 

If a paralysis of dynamics has occurred, a "mutation death" is evaluated as follows. TTte check 
for MRC hixel death relates to the verification of a suprathreshold resonance, where the 
30 resonance is defined as the integrand of partial resonances over the entire DST structure, ue. 
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Dsr 
>■ K 



DST DST, 



> 0 



DST 

0 



otherwise 



V r,.r,^ predef.ncdupperHmit.TypicallysettolOOforouteritcrationand 
1 000 for inner iterations. 

The time scale for realization of the Lindbald condition is changed and the system 
reiterated. 

H»oe fte final o«,«. of s.cp 314 is all hixcis ^ idantlly complmW 
o,ig„„udeoUde bWlngs in .he Mo^cl ,=»pl. which .-c ^ 
se,(h,(ij))whe«(h.W))istheoon«,»n<lin6olig=»»=leofti.se<r«nc t<»a,...ad for*. 

kth surviving hixel. 



p.cnn.nt Interactio" ^rientifv Diseases 

mutations identified using blodc 314 are processed using sim.lar steps w.thm block 
325 to identify diseasesrepresented by themutations-Hencethefmai output ofstep 335 .sas^ 

fLpressed!enestHatareassoc.ted.ithap^^^^^ 
.5 comJutatio„a.lybytheset(V.(y))v.here(V^^^^^^ [ooa-.a.] for the .-th element of the 

pathway capturing the k-th disease. 

1 • %^ ie block 335 can be omitted. 

For single disease analysis steps 324 - 334, i.e., diqck. o 
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riii'!trr' "P Prop ert y Check 
Di..«;id.n.inLi.i.lo*3;5 a,e p^cs^ed a, S..P 34. .0 id».ifr « P-"^ « 

iws. c,»«h,. "p-"- ^ — ' °"r 2 

fln^rtngofsuW^shold resonance and »bse^n»ec™te.„g(i.=.-»™va« of d.sp»s,^^^ 
Tho olns^Hng con.p„-on is hased on o,d»i„g (is h.s«. on „ 

o. . oiigo„n«d. -----r.TvpZ:^:;""-^- 

transversal of a hypergraph H - {X. E,, bz, tin.; 

T r^ E, .^forI=l,2,...,m,whereE,.E....EMdefmesubgraphs. 

1„ this method, each oligonucleotide, associated with a mutation that survives 'W- 
death" during resonant reverberation iteraUons. is represented by [aoa,...a.]. where a = 
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lA C TG) base associated is treated as a subgraph of the tota. set of unkrtown mutations that 
a,.' a'ct^ally present in the target sample. If the surviving hixel Is an ensemble than «.ch 
ensemble is treated as a subgraph with multiple nodes and several edges. If only an tndmdual 
hixel survives than it is treated as a single node subgraph. Transversal number of a hypergraph. 
H, is defined as the minimum number of vertices in a transversal. It is given by : 

f{H) = min^ 

The sub-steps involved in clustering are: 

determine MinD=U,^....,4}. -here A,. A. ... A. denote the sunriving 

resonance clusters. 

Next determine the following families: 

33= !3M^) ^'-^=^3^= Mm(rr3, vrrU}) 

3,=.... 

If Min A has k members, then the algorithm constructs Tr A = Tr in k steps. 
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Hashing Projector . u v tu^ 

A hashing p^l^ ,s .pplW « s»p 3« » «« ou-pu. of « check. Th 

hashing p^iec^r ^ » — i«» " "* "'f 

,0 h«c h.gh«. ^ '^■^ ^y^- """"^ " «" «>»^ 

distance of, at most two, are chosen, 
niao nnstic Decision 

A di^noalc decision is »tod a. s>cp 344 based upon d,c ontpu. of *e hash,ng ^y^. 
25 Th.drBnos,i.dcciston,s.chieved«singasi.p,..ab..ioo.upd»..indcxcdbyd„»sul>sof 

tashtaj ptojecaon computation using the aforementioned tables. 



Alternatives 

Alternative possible mutations are evaluated at step 346. 
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,f .te»dve, 3. available, *e a,W».iv. s« of n>~ions of i««s. a,. an- <^ 

"° . . . ,M H««« ittheoriginal«ofmutalk)iisfrom«hichtteonsiiial 

.„p«,edb=g™,«ga.sup3...H««^.««»^ 

,e.„fQ,F-sw..S=aen»d »n.,,^off^^ P ^ ^ ^^^^ 
alternative mutations, then the off-hne process is repe 

' ^"l'««<"«^o.y«dsMUo.^do.,--^.----^— ^ 
all rossible hypotees «. pcovkW as plausibk eandidatts. 

expressed genes and diseases. 

,5 co^ U.S. Pa«n. App..«ion 09/253.791. ».Ma- Txpon^-Ualb- Conv«^ 1^ 
M0.U0H.. V.1 Ua. M~»». ^ * l.=o,po™«a b, 



herein. 



■ exen-pla-y ««boa in,e„„ hav. been pmnaril, d«crib«. r.ftr«>c «. «ow Cham 

,l,.sJ^.pZ.l-«f.^-^"-^--*°- 
20 oft»roo,^n. fo, perform, ^ co«apo»dlns «ap. con^- 

::::r:::rpi: ..p jl™ of a ,«c« s^ 

ZL n, «ai.. R=*=r. * *os. compo.«.« n^ssary for a «,o,o„,h ™,de,s.a«d» of 
^!:llba«ba»i,l«s««aan.d«=Hb.dl.d«.».A-i.P.c— 
25 mo« compon»»s o,, «P0" «« i«pl«»cn«ion. .nay co„«in ftwe, "«"P«'"=»'»- 

^«.a„„^Ls.««p«s«— V.rtous..dl0..ions,o*.ee™^,™^^^^^^ 

^ ..dlly appa»«» s«M in *e ■« and .»eHc pnneiples defined h»e.n ma, ^ 
^i:^!ll=n.b«l.ne.«»i*o«,benseofUKinv.n.i.c.c«,.,T.^-— ^ 
3, ld«.»b..i»i«d»*ee.bo<.i™en.sh„w„he„inbu.is»b.a«^ed*ew.dcs.sco^ 
consistent «itl, Ihe principles and novel features disclosed he»m. 
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